Signatures of quantum stability in a classically chaotic system 
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We experimentally and numerically investigate the quantum accelerator mode dynamics of an atom optical 
realization of the quantum 5-kicked accelerator, whose classical dynamics are chaotic. Using a Ramsey-type 
experiment, we observe interference, demonstrating that quantum accelerator modes are formed coherently. We 
construct a link between the behavior of the evolution's fidelity and the phase space structure of a recently 
proposed pseudoclassical map, and thus account for the observed interference visibilities. 
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The relationship between the behaviour of classical and 
quantum systems, and how macroscopic classical phenomena 
originate in the quantum regime, remain subjects of dispute 
[fit]. The issues involved are particularly marked for quantum 
versions of classically chaotic systems [||. Experimental in- 
vestigations of such systems began with studies of microwave- 
driven hydrogen [Bfl; subsequent work has also centred on mi- 
crowave cavities H], mesoscopic solid-state systems and 
atom optics [^|, the approach we adopt. In this Letter we con- 
sider the quantum ^-kicked accelerator [|7l [| 0], a <5-kicked 
rotor with an additional static linear potential. The <5-kicked 
rotor is one of the most extensively investigated systems in 
chaotic dynamics [|lO||, and is equivalent to a free particle sub- 
jected periodically to instantaneous momentum kicks from 
a sinusoidal potential. Quantum mechanically, the effect of 
these kicks is to diffract the particles' constituent de Broglie 
waves into a series of discrete momentum states. In the 6- 
kicked accelerator, the linear potential modifies the chaotic 
classical dynamics only slightly, yet can radically change the 
quantum behaviour. The phases accumulated between con- 
secutive kicks by the momentum states are altered, leading to 
the creation of quantum accelerator modes [0, [| tm. We re- 
alize quantum 5-kicked accelerator dynamics in laser-cooled 
cesium atoms by the application of short pulses of a vertical 
standing wave of off-resonant laser light, which constitutes 
a sinusoidal potential; gravity provides the linear potential. 
Quantum accelerator modes, absent in the analogous classi- 
cal dynamics, are observed and are characterized by a linear 
(with kick number) momentum transfer to a substantial frac- 
tion (~ 20%) of the atoms. If coherent, this efficient mo- 
mentum transfer promises applications in atom interferometry 
Jl"i| ] . In this Letter we use a Ramsey-type interference exper- 
iment [|l2| to show that quantum accelerator modes do pre- 
serve coherence. We then relate the Ramsey fringe contrast 
to the fidelity / jl3[]; by a numerical analysis, we link the 
behaviour of / to the phase space structure of the <5-kicked ac- 
celerator in a pseudoclassical limit recently proposed by Fish- 
man et al. fll4{], Finally we explain differences in the observed 
fringe visibilities by examining the effect of the experimental 
range of kicking strengths. 

In our interference experiment, the atoms undergo <5-kicked 
accelerator dynamics, between the application of two tt/2 mi- 
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FIG. 1: Experimentally measured momentum distributions as the mi- 
crowave phase difference 9 is varied in a 7r/2 — 20 kick — 7r/2 
sequence where S b L ~ 35 GHz, with (a) T — 60.5/is (accelerator 
mode at -17HG), and (b) T = 74.5/xs (accelerator mode at 20HG). 
The corresponding numerically generated momentum distributions, 
where S b L = 35 GHz, are shown in (c) and (d). Population has been 
arbitrarily normalized to maximum value = 1. 



crowave pulses that couple two atomic hyperfine levels. In 
the absence of coherence-destroying spontaneous emission, 
the contrast of any interference fringes is related to the over- 
lap of two initially identical motional states evolved under the 
influence of slightly different Hamiltonians [jl5|], i.e. the fi- 
delity. It can therefore yield information on the sensitivity 
of the atoms' evolution to variations in the kicking strength. 
Strong sensitivity can be considered a quantum signature of 
chaos, particularly in the semiclassical limit (ft — > 0), hence 
the use by Peres JTrj ] of / as a measure of quantum stability. 

After magneto-optic trapping and molasses cooling to 5/iK, 
we prepare around 10 6 freely falling cesium atoms in the 
F = 3, nip = hyperfine level (denoted \a)) of the 6 2 Si/2 
ground state Jl7|]. The first n/2 microwave pulse creates an 
equal superposition of the atoms' internal states, i.e., \a) — > 
(|o) - ie ie \b))/V2, where \b) denotes the F = 4,m F = 
level. The phase 8 of this pulse can be changed with respect 
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to that of the second 7r/2 pulse, applied following 20 equally 
spaced 500 ns pulses from a standing wave of light. This 
is formed by retrorefiection of a Ti: sapphire laser beam; its 
maximum intensity is ~ 1 x 10 8 mW/cm 2 [[)]], and the light 
is red-detuned by 45 and 35 GHz from the Dl transition for 
atoms in states \a) and \b), respectively. After the second tt/2 
microwave pulse, we measure the momentum distribution in 
state | b) by a time-of-flight method. For more details of our 
experimental setup see Refs. [|| ^j. Measurement of a pe- 
riodic variation with 9 in the accelerator mode population in 
state \b), i.e. interference, directly implies coherent evolution. 
In the limit of large detuning, the Hamiltonian is 

H = H a \a)(a\ + H b \b)(b\ + ^(\b)(b\ - |o)(o|), (1) 
where %u> ab is the energy gap between \a) and \b), and 
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^[l + cos(GS)]5^<y(t-nT) (2) 



is the quantum ^-kicked accelerator Hamiltonian, acting on 
atoms in internal state \a) G {|a), \b)}. Here z is the ver- 
tical position, p the z-direction momentum, m the particle 
mass, g the gravitational acceleration, t the time, T the puls- 
ing period, il the Rabi frequency, t p the pulse duration, <5£ 
the detuning from the Dl transition for the state |<j), and 
G = 47r/ A, where A = 894 nm is the laser wavelength, and 
hG is a grating recoil (the momentum separation of adjacent 
diffracted states). We denote the amplitude of the phase mod- 
ulation to atoms in state \b) that results from application of 
the standing wave as (fid = Vl 2 t p /85 h L . The experimental 
mean value of (fid is 0.87T, and, due to the different detun- 
ing, that of the corresponding quantity for atoms in state \a) is 
'fid = 4>dh\l&i\ = 0.67T. We thus have effectively two differ- 
ent Hamiltonians, applied to the same initial motional state. 
The pulse train leads to the creation of a quantum accelerator 
mode, the momentum of which is the same for the two internal 
states [|[. We consider pulse periods T = 60.5/US and 74.5/zs, 
close to T1/2 — 2irm/frG 2 = 66.7/is, which corresponds to 
the lowest second-order quantum resonance in the <5-kicked 
rotor [|[ 11^]. For these T, well-populated accelerator modes 
involving substantial momentum transfer are created [0 H §]■ 
Figure 0(a) shows the measured final momentum distribu- 
tions of \b) atoms, for T = 60.5/^s. We see a period-27r vari- 
ation with 9 in the accelerator mode population (at around 
-17HG), the visibility V of which is (21 ± 2)% [Jo|]. We 
observe similar fringes for a range of detunings (8 L = 20- 
40 GHz) and total number of kicks N = 10-30. However, 
depending on their exact values, V can vary between 10% and 
40%. This periodic variation of the population demonstrates 
interference, and hence that the accelerator mode transfers 
momentum coherently. At T = 74.5/zs [Fig. 0(b)], however, 
fringes in the accelerator mode (at around 20hG) are prac- 
tically invisible, despite the expected coherent nature of the 
momentum transfer. In Figs. |l|(c) and 0(d) diffraction-based 
numerical simulations [0 ^ Q]7 incorporating the experimen- 
tal range of (fid (0.37T to 1.27r), also show this difference in the 



fringe visibility for the two values of T. The range of (fid is due 
to the Gaussian profile of the standing wave intensity (FWHM 
~ 1 mm) and the spatial extent of the atomic cloud (Gaussian 
density distribution, FWHM ~ 1 mm) [0] . As we optimized 
the overlap of the laser beams with the atomic cloud, the in- 
tensity and density maxima can be assumed to be coincident. 
The calculated visibility is then 25% for T = 60.5/is [Fig. 
|IJCc)] but only 8% for T = 74.5^s [Fig. 0(d)]. 

In order to explain these surprising observations, we intro- 
duce the Floquet operator Ff,((pd)- This describes the effect 
of one kick and the subsequent free evolution on the motional 
state of atoms in state | b) , where 

F b {(j) d ) = exp(~i[jx + p 2 /2]/k)cxp(i(j) d [l + cosx})- (3) 

We define F a (<f>d) analogously for state \a), with <fid replaced 
by (fid fll^]. As in Ref. [^], we use scaled position and momen- 
tum variables x = Gz and p — GTp/m, while 7 = gGT 2 
describes the effect of gravity, and k = TiG 2 T/m — —i[x, p] 
is an effective scaled Planck constant. After N pulses an ini- 
tial plane wave \q) of wavenumbergevolvestoFo-^d)-^^) = 
e t<pN^q (<^> d )) 5 where <fi — <fid or <fi a d , as appropriate. Regard- 
ing the initial motional state as an incoherent superposition of 
I q), the momentum distribution in state |6) for a given (fid after 
the 7r/2 — N kick — tt/2 sequence, is 

n(fep) = \ [dqC(q)[\r a (<t>d,p)\ 2 + \ip q b (<t>d,P)\ 2 ] 



dqC(q)r a (<f>d,pyTpt(<f>d,p) 



x cos(^/(^ d ,p) + NScfid + 0), 
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where 5(fi d = <fid - <fi a d , ifi q a {<fid,p) = {p\ip%(<fid)) , and </>/ is 
the phase of the interference term, i.e., j dqC (q)ifi1* ifil — 
I / dqC{q)ifi a *ifil\e % ' 1 " . The weighting C(q) describes the ini- 
tial Gaussian momentum distribution (FWHM = 6 grating 
recoils hG). The third (interference) term in Eq. (Q) is re- 
sponsible for the appearance of fringes in the accelerated \b) 
population. We denote the amplitude of the modulation in P b 
by A{dp d ,p)/2 = \$dqC{q)r a {4>d,p)*i>l{<t>d,p)\/% where 
J dpA(<fid,p) 2 — f(<fid) is the fidelity for a given (fid. 

We have calculated the individual terms of P b numerically 
for a wide range of (fid, obtaining as a consequence an im- 
portant result linking the pseudoclassical analysis of Fishman 
et al. Jl4] ] with the quantum stability measure of Peres JT^]. 
Comparison of Figs. ||(a) and |2|(c) with Fig. shows that the 
region in momentum space corresponding to an accelerator 
mode is also a region of high A. This remains high up to large 
values of (fid [pH, continuing beyond the point at which it has 
decayed to nearly zero in other regions of momentum space. 
As / = J dp A 2 , its large value when determined by integrat- 
ing over the momenta populated by atoms in the quantum ac- 
celerator mode implies that these atoms inhabit a stable region 
of quantum state space. Note that small A does not necessar- 
ily imply low atomic population, as can be seen in the plots of 
the non-interfering population J dqC(q)(\t/j^\ 2 + |?/^| 2 )/4in 
Figs. |(b) and |(d). Contrasting Fig. |(a) with Fig. |(c), we 



FIG. 2: Plots of A/2 against <j> d with N = 20, for (a) T = 60.5^s 
and (c) T = 74.5/is. Plots of the non-interfering population 
/ dqC(q)(\j>l\ 2 +\i>%\ 2 )/4 for (b) T = 60.5^s and (d) T = 74.5^s. 
Dashes mark the boundaries of the experimental range. Accelerator 
modes exist at — YThG for (a) and (c), and at 20HG for (b) and (d). 



see that this large value of A extends over a significantly wider 
range of (fid for T = 60.5/is than for T = 74.5/is. Hence, we 
can interpret the accelerator mode at T — 60.5/is as being 
more robust to variations in (fid, i.e. more stable, compared 
with that at T = 74.5 /is. However given our comparatively 
narrow experimental range of (fid, this does not explain the 
difference in fringe visibilities seen in Fig. |J 

The appearance of quantum accelerator modes in the 6- 
kicked accelerator is explained in the analysis of Fishman et 
al. Jl4| ] in terms of islands of stability centred on stable fixed 
points in the phase space generated by the map [|22|] : 
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Pn+i = Pn - ksm(xn) - sign(e)7, (5) 
Xn+i = X„ + sign(e)/3„ + i, (6) 

where the population of a mode is proportional to the size 
of the corresponding island. This is a pseudoclassical [e = 
(k — 2tt) — ► 0] rather than semiclassical (k — > 0) limit of 
the quantum dynamics characterized by the Floquet opera- 
tor of Eq. (^J). We have introduced p = pe/k (in an ac- 
celerating frame Jl4]]) and k — <fid\e\- Classically, the sys- 
tem is globally chaotic for our parameter regime. Figure || 
shows the pseudoclassical phase spaces generated by itera- 
tion of Eqs. (|^) and (^) for the experimentally investigated 
values of e = 2ir(T /T\m — 1), and a range of (fid- When 
(fid = 0.37T [Figs. ||(a) and ||(d)], the island (based around a 
stable fixed point in phase space) is substantially smaller for 
T = 74.5yus than for T = 60.5/is. For the average exper- 
imental value of (fid (0.87r) [Figs. ||(b) and||(e)], the islands 
have both grown to be about the same size. For (fid = 1.5ir 
[Figs. |3|(c) and 0(f)], the island has shrunk dramatically in the 
case of T — 74.5/is, while at T = 60.5/is the island has 
also shrunk, but not to the same extent. We therefore con- 
clude that the stable island representing the accelerator mode 
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FIG. 3: Stroboscopic pseudoclassical Poincare sections determined 
by Eqs. (|) and ^ for T = 60.5^s (=> e = -0.58), and for (a) 
(fid = 0.357T, (b) the average experimental value 0.87T, and (c) 1.57T. 
Figures (d), (e), and (f) show corresponding plots for T = 74.5/tts 
(=> e = 0.73). Displayed units are dimensionless. 



is much more robust to perturbations in the kicking strength 
for T = 60.5ns than for T = 74.5^s. The fact that A (and 
therefore /) remains large at the accelerator mode momentum 
for a significantly broader range of (fid when T — 60.5/is than 
for when T = 74.5/is, as shown in Fig. ||, matches exactly the 
observed greater stability of the island in the pseudoclassical 
phase space for T = 60.5/is. This is consistent with Peres's 
identification of the behaviour of the fidelity as reflecting sta- 
bility properties of the phase space in the semiclassical limit 
JTr^], even though our experiment (and numerics) operate in a 
pseudoclassical regime which is far from semiclassical. 

The position of the islands in pseudoclassical phase space in 
Fig. |3| indicates the region of the quantum accelerator modes' 
spatial localization. For T = 60.5p,s this is where there is zero 
standing light wave intensity, whereas when T — 74.5/is, it is 
where the intensity, and hence phase shift, are maximal [0]. 
We thus expect the modulation of the P\, interference term in 
Eq. (Q), cos((fii + N8(fid + 8), to have a strong dependence 
on (fid for the momenta at which accelerator modes are found 
when T = 74.5/xs, but not when T = 60.5/is. This is con- 
firmed by Figs. |(a) and |(b), where cos(0/ + N8(fi d + 8) 
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FIG. 4: Plots of cos(^j + NSfa) against <f> d with N = 20, for (a) 
T — 60.5/^s and (b) T — 74.5/is. Crosses mark the boundaries of 
the experimental range. 



is plotted as a function of p and 4>d for constant 9 (set to 
for convenience) for T — 60.5/is and T = 74.5/xs, re- 
spectively. At the accelerator mode momentum, the value of 
cos(4>i + N8<f>d) at T = 60.5/is is almost independent of <f>d, 
whereas at T = 74.5/is there is an approximate frequency 
doubling, relative to other momenta. We can now explain the 
presence or absence of interference fringes in Fig. [I] in terms 
of this difference in <fid dependence. As Eq. (|]) determines the 
population for one specific value of <f>d only, the observable 
momentum distribution is P&(p) = J dcf)dD{<t>d)Pb(<fid,p), 
where D((j>d) is a distribution describing the proportion of 
atoms experiencing a particular (f>d- Although for a single 
value of <j>d the visibility of the fringes at both T = 60.5/is and 
T = 74.5^s is high, Figs. |](a) and Q(b) show that integration 
over the experimental range of <pd causes a greater reduction 
in visibility at T — 74.5/is than at T = 60.5/is. At large <j>d 
(greater than occurs in our experiment), there is a breakdown 
of all structure in plots of cos(0/ + NStfid)- Comparing with 
Figs, ga) and g c), we observe this to coincide with a falloff 
in A. Note, however, that as <f>i is determined numerically 
from complex interference terms, one should be careful about 
attaching significance to values of cos(0j + N8<fid) for which 
the corresponding value of A is close to zero. 

In summary, we have performed a Ramsey-type interfer- 
ence experiment and thus demonstrated the coherence of the 
production of quantum accelerator modes, and hence their 
suitability for applications in atom interferometry. Numeri- 
cally, we have found the accelerator modes to correspond to 
regions of greater quantum stability, as quantified by the fi- 
delity. This is consistent with the presence of stable regions in 
the phase space of a recently proposed pseudoclassical limit 
of 5-kicked accelerator dynamics, rather than the globally 
chaotic behaviour of the semiclassical limit. These regions 
dictate the position of the accelerator modes' spatial localiza- 
tion, allowing us to explain the lack of fringes for the accel- 
erator mode at certain pulse periods, due to the experimental 
range of kicking strengths. Our investigation of coherence in 
quantum accelerator modes has allowed observation of their 
quantum-stable dynamics in this classically chaotic system. 
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